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ABSTRACT 

Particle probability hypothesis density filtering has become 
a promising means for multi-target tracking due to its ca¬ 
pability of handling an unknown and time-varying number 
of targets in non-linear non-Gaussian system. However, its 
computational complexity grows linearly with the number 
of measurements and particles assigned to each target, and 
this can be very time consuming especially when numerous 
targets and clutter exist in the surveillance region. Address¬ 
ing this issue, we present a distributed computation particle 
PHD filter for target tracking. Its framework consists of sev¬ 
eral local particle PHD filters at each processing element and 
a central unit. Each processing element takes responsibility 
for part particles but full measurements and provides local 
estimates; central unit controls particle exchange between 
processing elements and specifies a fusion rule to match and 
fuse the estimates from different local filters. The proposed 
framework is suitable for parallel implementation and main¬ 
tains the tracking accuracy. Simulations verify the proposed 
method can provide comparative accuracy as well as a sig¬ 
nificant speedup with the standard particle PHD filter. 
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1. INTRODUCTION 

Multi-target tracking(MTT) is to jointly estimate the num¬ 
ber of targets and position from a set of uncertain observa¬ 
tions. The classical approaches such as the nearest neigh- 
bor(NN) , Joint Probabilistic Data Association filter( JPDA) , 

and Multi-hypothesis tracking(MHT)[^ are based on the as¬ 
sociation algorithms. 
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Recently,much work has devoted to random finite set(RFS)- 
based approximations such as the probability hypothesis 
density (PHD) [3, [^,cardinalised PHD(CPHD)[ 8 ] [| and mul¬ 
tiple target multi-Bernoulli(MeMBer) filter]^. These meth¬ 
ods avoid the data association problem. From implemen¬ 
tation perspective, a particle PHD filter for nonlinear non- 
Gaussian MTT problems,was proposed in 0 And Vo 
proposed a close-form solution to PHD filter with assump¬ 
tions on linear Gaussian system. It is called GM-PHD filter. 

The demands of ’’real-time” MTT have been growing. Since 
the GPHD filter propagates both the intensity of the RFS 
and the posterior cardinality distribution,its real time char¬ 
acteristic is intrinsically not as good as the PHD filter.The 
MeMBer filter is more suitable for low clutter environments.And 
the GM-PHD is constrained to linear Gaussian system. The 
particle PHD filter is more suitable for nonlinear non-Gaussian 
MTT problems in dense clutter environment.However,the 
particle PHD filter computational complexity is very high,Therefore,we 
are interested in improving the real-time performance of the 
particle PHD filter. 

Similar to the particle filter,the resampling is chief obsta¬ 
cles to parallel and distribute. had proposed the dis¬ 
tributed particle resampling algorithm and implement in 
many WSN applications.In this paper,we proposed a dis¬ 
tributed particle PHD filter. In order to apply particle PHD 
filter theory to practice,distributed algorithm is better than 
centralized algorithm. Since the distributed algorithm can 
reduce the execution time by implementing the particle PHD 
filter using multiple processing elements. We distribute the 
particles to N PEs.Each PE can run particle PHD filter inde- 
pendently.All individual PEs compute local estimates based 
over observations in parallel,and transmit their estimate to 
a GU to obtain global estimate. 

Moreover,if the PEs were let to run as independent parti¬ 
cle PHD filter,each of them most likely have performance 
degradation caused by part of particles.The particles ex¬ 
change will solve the problem. 

The main contribution of this paper is summarized as fol¬ 
lows. First,we proposed the distributed particle PHD filter 
architecture. Second,we use the stphd method to extract 
the state estimation and corresponding measurement label. 

The real time performance is enhanced and the tracking 
performance is equal to traditional particle PHD filter. 


The remainder of this paper is organized as follows, stan¬ 
dard particle PHD filter is given in section 2.Section 3 ,we 
present our distributed particle PHD filter in detail.Simulation 
results are given in Section 4. Sect ion 5 is devoted to the con¬ 
clusions. 

2. BACKGROUD 
2.1 The PHD FILTER 

The PHD filter was developed in the framework of Finite 
Set Statistics(FISST) initially. The PHD function Ds is the 
first order moment of the random finite set(RFS) S and can 
be defined by 

Ds{x) = ElSsix)] = j 5x{x)Ps{dX) (1) 

where random density represen¬ 

tation of S. Ps is the probability distribution of the RFS 
S. The PHD has the properties that, the integral over a 
measurable subset S C E Ds{x)X{dx) is the expected 
number of target.In addition,the peaks of the PHD function 
give the estimates of the target states. 

PHD filter consists of the prediction step and the update 
step.Assuming the RFS is Poisson, it has been shown that 
the recursion propagating the PHD of fho multi-target 
posterior ph\k follows 

Dk\k — o ^k\k-i)Dk-i\k-i (2) 

where o represents composition of functions,4>fc|/e_i is the 
prediction operator and is the update operator,which are 
defined as follows: 


2) Update step 

For each z ^ Zk compute 

Lk-l + Jk 

Ck{z)= E C^) 

J = 1 

For i=l,...,Lfc_i + J/e update weights 

3) resampling step 

Compute the total target number Nk — resample 

{xi,wi/NktE''‘} to get 

Just like particle filter,the application of particle PHD 
filter is limited to its computational complexity which is 
mainly caused by the resampling and it also caused by the 
update which need all particles participate. 

2.3 Distributed particle filter 

The method of DRNA(distributed resampling with non¬ 
proportional allocation) was initially proposed by Bolic().The 
idea of the DRNA PF is to divide the whole particles into 
several groups so that the resampling can be performed inde¬ 
pendently by group and thus be implemented in parallel.The 
general DRNA is outlined by: 

1) Exchange particles among groups 

2) For k=l,...K and i=l,...,N sample x^ ~ 7r(xt) in par¬ 
allel in each group 

3) For k=l,...K and i=l,...,N compute the weights in each 
group in parallel 


($fc|fc-ia)(®) = 7fc + J A\k-i{x,Oo!{OHdO (3) 




1— pd{x) + 

zEZk 


i^k{z)E < ^k,z,o^ > 


for any function a on Eg,where 

0A:|fc — ^k\k — l{,^^ fk\k — l{,x\^^ T —1 

'ipk,z{x) = PD{x)gk{z\x) 
k{z) = XkCk{z) 

< f,g >=G f(x)glx)X(dx) 


a{x) 

(4) 


2.2 Particle PHD Filter 

As an approximate implementation of PHD filter,the par¬ 
ticle PHD filter is composed of three steps: 

At time k > 0,let Lk and Jk denote the number of survival 
particles and birth particles at time k,repectively 
1)Prediction step: 

For i=l,...,Lfc_i,sample x^. ~ qk{-\xk-i, Zk) and compute 
the predicted weights 


^k I A: — 1 


4 ^k\k — l{Xk^ Xk — i) I 

qk{xi\xl_^,Zk) 


(5) 


For i=L/e_i + l,...,Lfc-i + Jfc,sample xl ~ pk{-\Zk) and 
compute the weights of new-born particles 


^k I A: — 1 


Ihixl) 1 
Pk(yX^Zk^ Jk 


( 6 ) 


— - Tirx - (9) 

4) Normalize the weights of the particles with the sum of 
the weights in the group 

5) Resample inside the groups 

3. DISTRIBUTION PHD FILTER 

The probability hypothesis density(PHD) filter,which prop¬ 
agates only the first moment instead of the full multi-target 
posterior,still involves multiple integrals with no closed forms 
generally(by vol).So BA-NGU VO proposed the particle PHD 
filter.The particle PHD filter is suitable for problems that 
nonlinear non-Gaussian dynamics.However,it’s high compu¬ 
tational is to hold back it’s application into real time sys¬ 
tem. To speed up the particle PHD filter, we propose a 
distributed approach that only uses a subset of particles to 
different computing cores. In other words,we use a subset 
of particles to different PE and acquire the same accuracy 
as all the particles were used together ,but avoid unneces¬ 
sary communication among the PEs. We entitle this idea 
DGPEPHD in this paper,can be formalized as follows. The 
structure chart of DGPEPHD algorithm is shown in Eig.l 

As the particle PHD filter is similar the particle filter,also 
involve three basic steps,we can use DRNA in particle PHD 
filter too.In this paper,we apply the DRNA scheme to particle- 
PHD filter.However,there is some differences in how calcu¬ 
late the weight and estimate target,the challenge is change 
the DRNA so we can use it in particle PHD filter(Eig 1). 










Figure 1: the structure of DCPFPHD 


3.1 General structure 

The standard particle PHD filter requires all the particles 
be participated by a (single processor) in resampling step. 
In this paper,Assume we have K PEs and each can run a 
separate particle PHD with M particle.The total number of 
particles is N=MK. The distributed computation particle 
PHD filter can be outlines as follows, 

Assume we have K groups and each can run a separate 
particle PHD with M particle.The total number of particles 
is N=MK. The distributed computation particle PHD filter 
can be outlines as follows. 


Algorithm 1 Depiction on Each Group j 
DCPFPHD 

Input: ,Z{k) 

Local PF-PHD filter 
for z = 1,..., M do 

end for 

for z = M + l,...,M + J/c do 

^p{-\Zk). 

end for 

for z = 1,..., M + Jk do 
weight update 

4"'' = Hxi) + E.ez, 

end for 

Local Estimation 

estimate Nk{j) and state(j) 

Resample 

resample the M Jk particles into Nk{j)Rk particles and 
set M = Nk{j)Rk 

Share Particles 

each group send t particle-weight couples to neighbors 


In the following,we describe the algorithm steps in detail. 
t=0:Generate M particles for each group.The number of 
all particles is N.All these particles are shared with same 
weights ujQ — 1/M 
t=tT 1 

Step 1: Predict ion At time t-l,we assume the particle set 
is available.For k-th group and m=l,...Lfc,sample 
1 1 YoT new-born particles we divide these 


particles to K group equal then join the groups. 

Step 2:Update In update step,we can calculate the weights 
among groups.Let Zk denote the measurements set.For each 
group,for each zinZk^nse the likelihood and compute Cf and 
then update the weights. 

Since we use the STPHD method to estimate state,the 
update step should make some change. zhao[^ has demon¬ 
strated that the Dk\k can be calculated as 


Dk\k{x\Zi.,k) = ADk\k{x\z) + ADk\k{x\(p) (10) 

zEZk 

where ADk\k{x\(l)) denotes the PHD of the measure unde¬ 
tected. For each observation Zp G Zk = !,•••, Mfc,For each 
group j,j=l,...,K 

M+J/g 

Ck{zk,p) = ( 11 ) 

i=l 


l/k,z{Xk^) 


( 12 ) 


^('^(fc,^) T Cki^Zk^p) 
then calculates the sub-weight of each particles for all ob¬ 
servations Zk,p 


Awk = 


(13) 


And the particle sub-weight for the target without measure¬ 
ments obtained is 

V p — i,j\ — i,j 


(14) 


Based on formula ([T0| ,the weights can compute through 

(15) 




p=i 


i,p,0 

k 


Step 3.'Local estimation. 

We use the STPHD method can extract the estimate tar¬ 
gets and these targets’ sequence number. For each measure- 

- '^'^■>Pt 3 

ment Zk,p,p — !,•••, M^,compute the sum of Awk relevant 
to Zk,p in group j. 


_ __ , 

Y. A^T' (16) 

i=l 

Gompute the sum of sub-weight corresponding to 

targets without observations: 

M+Jk 

= Y ( 17 ) 

i=l 

Since the weight sum of all the particles equals to the 
target number,the local target number can be estimated by 
LT/ = Wk^) where Tound{^l~^ Wk^) is the 

nearest integer to X) 

Find the LT/ largest sum weight AW//^ and the index set 
relevant to AW//^ . The local estimated target state can 
be calculated by where I in index set , 




~i,h3 

Awk 




( 18 ) 

















The groups can send data to the CU. When the group get 
a local estimate like Ck,i and transmit the pair {Ck.h^i} to 
the CU.Then the CU can combine the local estimate state 
which depends on the measurement index from the group to 
construct a global estimate. 

Step 4: Global Estimation. 

As the CU receives all the local information -Depend 

on the rule ”at most one measurement per target” [TU. If the 
local estimate state’s label {/j} is same that from different 
groups,then we use their mean value as the global estimate 
state. And the local estimate state may be from clutter,so 
we consider only the local estimate states’ number which 
have same label greater than half of groups’ number as valid 
estimate states. 

Step 5.’Resampling. 

The resampling can be carry out locally at the N groups.Norms 
the weights of the particles with the sum of weights in the 
group. 

Step 6:Local exchange. 

The particles in the n-th group will degenerate when its 
aggregated weight becomes negligible relative to the aggre¬ 
gated weights of the other groups.Then the n-th group hardly 
contribute to the approximation of the posterior probability 
distribution.In order to keep the groups are valid,neighboring 
groups can exchange a portion of particles and weights.We 
select the L particles from k-th group(L < M/2) to replace 
the L particles from k-1 group in random, 
for k=l,...K-l,i=l....,L do: 


cording to the following model Xk+i = 


T 

1 

0 

0 


0 

0 

T 

1 


Xk~\- 


TV2 


0 

0 

TV2 

1 


Wk where Xk = [xk.Xk.yk.Vk] is target 


state vector at time kT(k is the time index and T=1 is 
the sampling period).[xk,yk] is the position,while [xk,yk] 
is the velocity.rcfc = is the vector of independent 

zero-mean Gaussian white noise with standard deviations of 
■ 0.025 0 ] _ . . ^ 

Q ^ ihere is just a signal sensor m the scenario 

and the target-originated measurement are given by 
Zk=g{xk)-\rVk=[R, 


R= \/{x- SxY + {y- SyY 

-U 

1 

0 = arctan {y — Sy)/{x — Sx) 


. '^k 


for k=K ,i=l....,L do: 


The measurement variance Vk = 5m, vl = 0.05rad. Clut¬ 
ter is modeled as a Poisson RFS Hik with intensity Kk{z) = 
AfcUn(^) The Target can disappear or appear in the scene at 
any time.The probability of target survival is ek\k-i — 0-9 
and is detected with probability Pd,^ = 1 . Assume the 
target birth according to the Poisson distribution with the 
intensity N(-,x,Q) where N(-,x,Q) denotes a normal den¬ 
sity with mean x and covariance Q. 

f 0 \ f 10 0 0 0 \ 

0 10 0 

0 0 10 0 

\ 0 0 0 1 y 

The surveillance region is [—7r/2, 7r/2] x [0, 200] rad.m. As¬ 
sign 200 particles to an exist or new born target in each 
group.The number of group is 4. 


3 
0 

-3 / 


Q = 


4.1 Simulation Results 


The processors are connected using an interconnection 
network.There are many type of network can be used,we se¬ 
lect a ring configuration in this paper.For a ring network ,the 
nth PE,n=l,...,N-l,can send data to the (n+l)th PE.The 
Nth PE transmits data to PE number 1. 

The PEs can also send data to the CU. When the PE 
get a local estimate like ddd and transmit the pair to the 
CU.Then the CU can combine the local estimate state which 
depends on the measurement index from the PE to construct 
a global estimate. 


3.2 Analysis of Time Delay and Computational 
Complexity 

In the traditional particle PHD filter,all the particles have 
to be involved by serial. In our methods,since particles are 
divided into K groups and the group can run a particle PHD 
filter independently,thus it only utilise 1/K time than before 
in theory. 


4. SIMULATION RESULTS AND EXPERI¬ 
MENTAL STUDY 

To evaluate the proposed distributed particle PHD fil¬ 
ter, we consider a two-dimensional scenario with the target 
can disappear and appear at anytime.Each target moves ac- 


We run 100 independent simulations of the HOT model 
given by....The overall number of particles was N = 2000 
and ,for the DCPEPHD,we divide them into N=4 proces¬ 
sors with M=500 particles each. The tr ue tr aje ctorie s of 
five tracks over 50 scans are plotted in Eig |2(a) Eig 2(a) [ also 
shows the positions of the estimated targets over 50 time 
steps. The individual x and y coordinates of the track s and 
estimat ed ta rgets for each time step are shown in Eig 2(b)| 
and Eig |2(c) [ respectively.lt can be seen that estimated posi¬ 
tion based on the traditional particle PHD filter which the 
number of particles is equals the DPHD filter’s particles and 
the DCPPHD are similar and they are all close to the true 
tracks. 

Another vital factor for the performance of the method is 
the estimated number of targets.The number of true targets 
and estimated target by our method ,standard particle PHD 
filter with same particles and only 1/k particles at each scan 
are given in Eigj^It is observed that,under the same simu¬ 
lation conditions and same particles ,the particle PHD filter 
and DCPPHD filter achieve 

It was proposed in to use the Optimal Sub-Pattern As¬ 
signment (OS PA) as a multi-target miss-distance metric,and 
the parameters in it are set as p = 1 and c = 100 inour eval¬ 
uation. Eig 1^ shows the OSPA distance of the DCPPHD 
filter and particle PHD filter. 

The total number of particles for these methods is given in 
Tabell.Implemented on a Dell computer using MATLAB,this 
approximate times used by these approachs are also given 




























Table 1: 12345 


method 

particle number 

times 

mean of OS PA 

std of OSPA 

r 

PHD 

1 in 1,000 

17.6540 

3.0577 

0.0548 

0 

DCPPHD 

1 in 5 

8.6331 

3.0732 

0.0830 

0 

partPHD 

1 in 40,000 

4.0935 

3.1993 

0.2590 

0 

PHD 

1 in 1,000 

72.2214 

6.4432 

19.0788 

10 

DCPPHD 

1 in 5 

27.1516 

6.2220 

5.4630 

10 

partPHD 

1 in 40,000 

17.2492 

7.4752 

54.4956 

10 

PHD 

1 in 1,000 

117.1504 

8.0049 

8.2928 

20 

DCPPHD 

1 in 5 

39.4615 

8.0884 

22.6434 

20 

partPHD 

1 in 40,000 

29.2382 

9.7143 

136.3478 

20 



(a) 



(b) 



(c) 


Figure 2: Target trajectories estimated by the DCP- 
PHD filter and the particle PHD filter.(a)Estimated 
trajectories and true trajectories, (b)Estimated tra¬ 
jectories and true trajectories in x-axis direc¬ 
tion. (c)Estimated trajectories and true trajectories 
in y-axis direction 
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Figure 3: Estimated number of targets and the true 
number of targets 
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Figure 5: the time 

in Tabell.And the Tablet shows the mean of OSPA as well. 

set simulation results l.the target position 

2. the ospa 

3. the time 

4. the number 

and analyse the result. 

5. CONCLUSION 

In this paper,we proposed a DRNA particle PHD filter 
that want to improve the particle PHD filter runtime.The 
DRNA-PHD filter can speed up the particle PHD filter in 
theory.However, the feasibility of the proposed method needs 
to be tested in real applications.lt note that divide the more 
groups will lead to the decrease of performance, 
future work may consider other particle exchange method 
and simplify the update step as we found the weight update 
also cause a lot of compute time.The update step and re- 
samping will be the bottleneck for the development of higher 
speed of particle PHD filter.And we will test the method into 
the GPU as well. 
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